Lattice study of the chiral magnetic effect in a chirally imbalanced matter 
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We investigate the chiral magnetic effect by lattice QCD with a chiral chemical potential. In 
a chirally imbalanced matter, we obtain a finite induced current along an external magnetic field. 
We analyze the dependence on the lattice spacing, the temperature, the spatial volume, and the 
fermion mass. The present result indicates that the continuum limit is important for the quantitative 
argument of the strength of the induced current. 



PACS numbers: ll.15.Ha, 12.38.Gc, 12.38.Mh 



O 

(N 
o 

CD 

Q 
o 



> 
00 



I. INTRODUCTION 



Topology plays a significant role in classical and quan- 
tum field theory. In quantum chromodynamics (QCD), 
the nontrivial topology of the background gauge field is 
related to the axial anomaly of the fermion. From the 
Atiyah-Singer index theorem. 



NfQ = Nr^Nl 



(1) 



Q. In the QCD vacuum, the total numbers of the left- 
handed and right-handed fermions are the same, and the 
average topological charge is zero. However, the topolog- 
ical charge strongly fluctuates in space-time. This topo- 
logical fluctuation is essential for the ry'-meson mass [1, Q . 

The chiral magnetic effect is one possible candidate to 
detect the topological fluctuation in experiments 0-0) ■ 
A noncentral heavy-ion collision produces a very strong 
magnetic field. The fermion flows along this magnetic 
field, and its direction is determined from its chirality. 
The topological fluctuation generates a chiral imbalance 
in some local domain, and then it induces a charged flow 
of the fermion, or equivalently, an electric current. When 
all experimental events are averaged over, the total elec- 
tric current is zero, but its event-by-event correlation is 
nonzero. 

Several theoretical works have studied the chiral mag- 
netic effect in lattice QCD 8-12j. These works have tried 
to measure current-current correlation or charge density 
distribution. In order to analyze the chiral magnetic ef- 
fect in this direction, it is necessary to take into account 
topological objects on the lattice, however, it is usually 
difficult in practice. We overcome the difficulty by intro- 
ducing chiral chemical potential, instead of the topolog- 
ical fiuctuation. 

In this paper, we investigate the chiral magnetic effect 
in lattice QCD by extending a previous work ^13J. We 
discuss a theoretical background of the chiral chemical 
potential in Sec. II, and the lattice QCD formalism in 
Sec. HI. We show the numerical result of the full QCD 
simulation in Sec. IV and of the quenched QCD simula- 
tion in Sec. V. Finally, we conclude in Sec. VI. 



II. CHIRAL CHEMICAL POTENTIAL 

The chiral chemical potential /^s is defined in the con- 
tinuum Dirac operator as 



Dills) = 7m (^M 



■ igf'Al) + m + ^57475 



(2) 



[14| . The Euclidean metric is used throughout this paper. 
The chiral chemical potential couples to the chiral charge 
density 
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= d^xiij^^R - V'I^l) = ^{Nr - Nl) .(3) 

The chiral charge is the number difference between the 
right-handed and left-handed fermions. Strictly speak- 
ing, the chiral charge is not a conserved quantity in QCD 
because of the axial anomaly. The chiral chemical poten- 
tial realizes a finite chiral charge in an equilibrium state, 
namely, a chirally imbalanced matter. This is useful for 
theoretical treatment, for example, the lattice QCD sim- 
ulation is possible. 

The chiral chemical potential has been introduced for 
analyzing the chiral magnetic effect [l3 - [l7t . The induced 
vector current was derived. 
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from the one-component Dirac equation coupled with the 
background magnetic field [3]. (Note that the overall 
factor is different from its original expression by q. This 
is a matter of definition. The electric current is usually 
defined as Jem = qJ- ) This formula is independent of 
the fermion mass and the temperature. 

The significant property of the chiral chemical poten- 
tial is that the sign problem does not arise. In lattice 
QCD, a quark chemical potential causes the sign prob- 
lem. The naive Monte Carlo algorithm breaks down at 
a finite quark number density. Although many numeri- 
cal methods have been proposed for simulating the quark 
chemical potential, their applicabilities are limited within 
a small quark chemical potential For a large chem- 
ical potential, we can only simulate exceptional cases. 



2 



such as the two-color QCD and an isospin chemical po- 
tential, which can exactly avoid the sign problem |19l - l23| . 
The chiral chemical potential is an interesting possibility 
to study the finite density matter in lattice QCD. 



TABLE I: Simulation parameters of full QCD ^2^. 



mps / my Ns Nt 



5.32144 1.5 GeV 0.16650 



0.5 
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III. LATTICE QCD FORMALISM 

For the numerical simulation, we used the SU(3) pla- 
quette gauge action and the Wilson fermion action. We 
considered the degenerate two-flavor case, in which the 
two fermions have the same mass m and charge q. This 
approximation simplifies the simulation algorithm, espe- 
cially, the hybrid Monte Carlo algorithm. 

The lattice Dirac operator of the Wilson fermion is 



i 

4(1 

(1 -I- 746 



[(l-707^.+ + (l+7.m-] 



746 ■ " )14+ 

1^575 "J 



with 



g±a/i575 ^ cosh(a^5) ± 75 sinh(a^5) 



(5) 

(6) 
(7) 
(8) 



This Dirac operator reproduces the continuum form ([2]), 
apart from an overall factor, in the naive continuum limit 
a — >■ 0. This Dirac operator satisfies the "75-Hermite" 
property 



or equivalently. 



(9) 



(10) 



From this property, the fermion determinant is real and 
semi-positive in the two-flavor case. 



det 



Difi5) 
D{fi5) 



detL>(^5)P > 



(11) 



Thus, the sign problem does not occur. 

To apply an external magnetic field, we also need 
the U(l) gauge field. On the lattice, the U(l) gauge 
field is given as the Abelian phase factor u^i{x) — 
exp(iaqAf^{x)) . Since the magnetic field is external, the 
field strength term of the U(l) gauge field is not intro- 
duced. The U(l) gauge field is introduced only in the 
Dirac operator ([5]) by replacing 



Uf_,{x) u^{x)U^{x) 



(12) 



In a finite-volume lattice with periodic boundary condi- 
tions, the quantized value of the magnetic field is allowed 

as 



a^qB 



27r 
iVF 



X (integer) 



(13) 



24| . For a homogeneous magnetic field in the 2:3- 
direction, the phase factors are 



ui{x) 

U2{x) 
Uf,{x) 



exp{—iaqB NsX2) at xi : 

exp(iagi?xi) 

1 for other components 



aiV, 



(14) 
(15) 
(16) 



The magnetic field B and the charge q do not appear 
separately. The combination oF'qB is an input parameter 
of the simulation. 



IV. FULL QCD ANALYSIS 

In this section, we show the numerical results in full 
QCD. To generate the dynamical gauge ensembles, we 
used the hybrid Monte Carlo algorithm. The lattice spac- 
ing a^^ ~ 1.5 GeV (a ~ 0.13 fm) and the mass ratio of 
pion to /9- meson is mpg/mv — 0.5 (25| . The simulation 
parameters are listed in Table ID The physical tempera- 
ture is T = l/{Nta) ~ 400 MeV, which is in the decon- 
finement phase. 

In Fig. [U we plot the chiral charge density 



"5 



(17) 



site 



At a very large chiral chemical potential, the so-called 
saturation, which is a lattice artifact, appears [13] ■ To 
avoid this problem, we restricted the chiral chemical po- 
tential in a/X5 < 0.5. A finite chiral chemical potential 
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FIG. 1: The chiral charge density ns in full QCD. 
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FIG. 2: The vector current density J as a function of the 
magnetic field qB in full QCD. 
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TABLE II: Simulation parameters of quenched QCD [26 



p 






mps /mv 


Ns 


Nt 


5.90 


1.9 GeV 


0.15920 


0.4 


12, 16, 20 


4 


5.90 


1.9 GeV 


0.15890 


0.5 


12 


4 


5.90 


1.9 GeV 


0.15740 


0.7 


12 


4, 6, 12 


6.25 


3.1 GeV 


0.15115 


0.7 


18 


6, 10, 18 


6.47 


4.0 GeV 


0.14885 


0.7 


24 
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density as 



(19) 



The factor A^dof = Nc x Nf = 6 is the number of 
fermions with the same charge. This functional form 
is consistent with the analytic formula ([4]). In the an- 
alytic formula, the overall coefficient is l/(27r^) ~ 0.05. 
On the other hand, the best-fit value of the lattice data 
is C = 0.013 ± 0.001. The induced current seems to 
be suppressed. However, in order to compare the over- 
all coefficients quantitatively, we have to estimate several 
systematic effects, e.g., the renormalization. In the next 
section, we analyze such systematic effects in quenched 
QCD. 



V. QUENCHED QCD ANALYSIS 

We performed the quenched QCD simulation. The 
quenched approximation is to neglect the fermion deter- 
minant in the Monte Carlo sampling. In quenched QCD, 
the computational cost is highly reduced and the param- 
eter tuning is very easy. This approximation is consid- 
erably reasonable in many cases, e.g., hadron spectrum. 



FIG. 3: The vector current density J as a function of the 
chiral chemical potential ^5 in full QCD. 



generates a finite chiral charge density. The chiral charge 
density is enhanced when the magnetic field is applied. 
This is because the thermodynamic potential is increased 
by the magnetic field [T^ |. 

We calculated the local vector current density 



J 
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V 



(18) 



along the magnetic field. The transverse components are 
always zero, (ipjiip) = {'4'l2'4') — 0, because they are 
irrelevant for the chiral magnetic effect The vec- 

tor current density is plotted as a function of the mag- 
netic field in Fig. [2] and of the chiral chemical potential 
in Fig.O The vector current density is a linearly increas- 
ing function both of the magnetic field and of the chiral 
chemical potential. We parametrize the induced current 




FIG. 4: The chiral charge density ns in quenched QCD. The 
physical temperature is about T ~ 500 MeV. The meson mass 
ratio is mps/mv — 0.4. 
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FIG. 5: The vector current density J in quenched QCD. The 
simulation parameters are the same as in Fig. [4] 



FIG. 6: Lattice spacing a dependence. The physical tem- 
perature is about T ~ 500 MeV. The meson mass ratio is 
mps/mv ^ 0.7. 



We expect the quenched approximation to work at least 
for a small chiral chemical potential, and evaluate sys- 
tematic effects on the induced current. The used gauge 
configurations are summarized in Table |lT] [26] . 

We show the chiral charge density in Fig. |4] and the 
vector current density in Fig. [5l The qualitative behavior 
is consistent with the full QCD simulation. The induced 
current shows the linearly rising behavior. We fit the 
induced current by Eq. ()19p and analyze the parameter 
dependence of the overall coefficient C. 

First, we focus on the problem of the discretization. 
The local vector current is not renormalization-group in- 
variant on the lattice. This is different from the contin- 
uum theory, in which the local vector current is strictly 
renormalization-group invariant. The conserved vector 
current on the lattice is a point-split-form Green func- 
tion [27| . To obtain the renormalization-group invari- 
ant forrn, we have to calculate the renormalization fac- 
tor [H, [2§|- In general, a large statistics is needed to 
calculate the renormalization factor of the fiavor-singlet 
vector current because the disconnected contribution is 
rather noisy. Instead of this, we analyze the dependence 
of the overall coefficient on the lattice spacing a. In the 
continuum limit a ^ 0, the renormalization factor au- 
tomatically approaches to 1, and thus the ambiguity of 
the renormalization factor disappears. The lattice spac- 
ing dependence is also important to estimate other lat- 
tice discretization artifacts. For example, the 0(a) term 
of the Wilson fermion action explicitly breaks the chiral 
symmetry. Such lattice artifacts also disappear in the 
continuum limit. 

We calculated the induced current at /3 = 5.90, 6.25, 
and 6.47. By tuning the simulation parameters, we 
fixed the spatial volume, the temperature, and the me- 
son masses. In Fig. [6l we plot the best-fit value of the 
overall coefficient C. The overall coefficient depends on 
the lattice spacing. Its value is close to the full QCD 



O 



0.05 



0.04 



0.03 



0.02 - 



0.01 - 



T r 



P = 5.90 I— e- 
P = 6.25 I— B- 



100 200 300 400 500 600 
T [MeV] 

FIG. 7: Temperature T dependence. The phase transition 
temperature is about 270 MeV. The meson mass ratio is 
mps/mv ^ 0.7. 



result in a > 0.1 fm and larger in a < 0.1 fm. This result 
indicates that it is important to reduce the discretiza- 
tion effect and to take the continuum extrapolation for 
more quantitative argument. Although we cannot deter- 
mine the functional form of the extrapolation, the overall 
coefficient seems C ~ 0.02 - 0.03 in the continuum limit. 

In Fig. [71 we plot the data of /3 = 5.90 and 6.25 as 
a function of the temperature T — l/{Nta). Although 
the unrenormalized values need not be the same between 
the different lattice spacings, the qualitative behaviors 
should be consistent. The data of /3 = 5.90 and iVj = 4 
(T ~ 500 MeV) deviates from the others. The reason 
would be that the discretization artifact of the unim- 
proved fermion is sizable at Nt = 4, which is known in 
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FIG. 8: Volume V dependence. The physical temperature is 
about T ~ 500 MeV. The meson mass ratio is mps/mv — 0.4. 



FIG. 9: Hopping parameter k dependence. The dotted line 
corresponds to the critical hopping parameter Kc = 0.15983. 
The physical temperature is about T ~ 500 MeV. 



the standard lattice QCD at finite temperature. Except 
for A^t = 4, the induced current is not so sensitive to the 
temperature. The induced current is nonzero even below 
the phase transition temperature Tc ~ 270 MeV, i.e., 
in the confinement phase. Basically, the chiral magnetic 
effect is expected in the deconfinement phase because col- 
ored particles cannot flow independently in the confine- 
ment phase. Such a real-time information is, however, 
not reflected in the local value of the induced current in 
an equilibrium state. 

Next, we show the volume dependence in Fig. El The 
lattice simulation is performed in a finite box, not in 
the infinite volume. The spatial volume should be large 
enough that the obtained result is independent of the 
spatial volume. In order to change the spatial volume 
V = a^N^, we fixed the lattice spacing a and changed 
the spatial extent as iV, = 12, 16, and 20. As expected, 
the induced current is independent of the spatial volume. 

Finally, we analyze the dependence on the hopping pa- 
rameter, which is equivalent to the fermion mass depen- 
dence. The result is shown in Fig. [HI The used hopping 
parameters are k = 0.15920, 0.15890, and 0.15740, which 
correspond to mps/my = 0.4, 0.5, and 0.7, respectively. 
The critical hopping parameter is Kc = 0.15983, which 
corresponds the chiral limit [1^. The induced current is 
almost independent of the fermion mass. This is consis- 
tent with the analytic formula @. 



VI. SUMMARY AND OUTLOOK 



quantitative level, the induced current is somehow sup- 
pressed compared to the analytic formula. The estimated 
value of the overall coefficient is C ~ 0.02 - 0.03 even af- 
ter the continuum extrapolation. Although the system- 
atic effects have been estimated in quenched QCD, the 
situation will be more or less the same in full QCD. In 
order to reduce the discretization effect, we should first 
introduce the improved fermion action, since the fermion 
action has a larger discretization error than the gauge 
action. More quantitative calculation will shed light on 
possible QCD corrections to the chiral magnetic effect 

We have not discussed about chiral symmetry. Since 
the Wilson fermion explicitly breaks chiral symmetry at a 
finite lattice spacing, we cannot analyze chiral symmetry 
in the present calculation. We show the functional forms 
of other lattice Dirac operators in Appendix A. By using 
lattice Dirac operators with better chiral property, we can 
study the role of chiral symmetry for the chiral magnetic 
effect. 
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We have analyzed the chiral magnetic effect by the lat- 
tice QCD simulation with the chiral chemical potential. 
At the qualitative level, the induced current is consistent 
with the analytic formula which is derived from the Dirac 
equation coupled with an external magnetic field. At the 



Appendix A: Lattice Dirac operators 

The simplest choice of the lattice Dirac operator with 
the chiral chemical potential is the Wilson Dirac operator 
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([5]). The chiral chemical potential can be implemented 
for other lattice fermions with better chirality, e.g., the 
staggered fermion, the domain- wall fermion and the over- 
lap fermion. 

In the staggered fermion, the Dirac spinor is con- 
structed by mixing spinorless Grassmann fields on differ- 
ent lattice sites. For example, in four dimensions, a four- 
taste four-component Dirac spinor is constructed from 
2"* lattice sites. The gamma matrix appears as a direct 
product (ja^lb) of two matrices which act on the spinor 
space and the taste space, respectively. A naive choice of 
the staggered Dirac operator is 

-Dks(M5) = ma+^^ ili{Ti+ - Ti_) 



+ -ry4(T4+e'"^^i^-^ -r4-e^"^^^^) , (Al) 



with the staggered phase factor 



(A2) 



We consider two types of T^. To construct the taste- 
singlet matrix (75 ^1), we take 



(A3) 



sym 



where "sym" means the symmetric average of the path- 
ordered products. The matrix factor e±°M5r5 jg ^ product 
of hopping terms, and thus nonlocal. This is understood 
from the fact that the operation of the gamma matrix cor- 
responds to mixing different lattice sites in the staggered 
fermion. Such a nonlocal Dirac operator is expensive for 
the numerical simulation, especially for the dynamical 



simulation. For a local realization, we take 

[T^Uy = [V5h.y = (-1)^1+^^+^^+^^,,, . (A4) 

This matrix is converted into the taste-nonsinglet ma- 
trix (75 ® 7^). The generator (75 (g) 7J) defines the 
U(l)xU(l) residual chiral symmetry, i.e., the exact chiral 
symmetry of the massless staggered fermion. The above 
formulation is based on four tastes, i.e., four flavors. To 
implement it in the 2-flavor or (2 -1-1) -flavor dynamical 
simulation, we need to take square root or fourth root of 
the Dirac operator. The subtle problem arises especially 
in taking the fourth root. Moreover, since this chiral 
chemical potential is not taste-singlet, the interpretation 
of the fourth root is nontrivial. 

Since the Wilson Dirac operator ([S]) is 75-Hermitian, it 
is straightforward to formulate the domain- wall fermion 
and the overlap fermion. This is different from the case of 
the quark chemical potential, in which the sign function 
of the non-Hermitian matrix must be introduced. The 
domain- wall Dirac operator is 



-Ddw(M5 



£'w(M5) 
1-7E 



+ 1 

^ 1 + 75 ^ 
-'5+ 7i J5- 



(A5) 

where is similarly defined to Eqs. © and (O in the 
fifth dimension but without the gauge field. The overlap 
Dirac operator is 



DovifJ-5) = 1 +75e(75^w(M5)) , 



(A6) 



where e{H) is the sign function of a Hermitian matrix H. 
These Dirac operators hold the expected properties, such 
as the 75-Hermite property and the Ginsparg- Wilson re- 
lation. 
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